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Abstract 

The interaction between intersubband resonances (ISRs) and metamaterial microcavities can form 
a strongly coupled system where new resonances form that depend on the coupling strength. Here 
we present experimental evidence of strong coupling between the cavity resonance of a THz meta¬ 
material and the ISR in a high electron mobility transistor structure with a triangular confinement. 
The device is electrically switched from an uncoupled to a strongly coupled regime by tuning the 
ISR with epitaxially grown transparent gates. An asymmetric triangular potential in the het¬ 
erostructure enables ultrawide electrical tuning of ISR which is an order of magnitude higher as 
compared to the equivalent square well. For a single triangular well, we achieve a coupling strength 
of 0.52 THz, with a normalized coupling ratio of 0.26. 
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Ultrastrong light-matter interaction is one of the key aspects of cavity quantum elec¬ 
trodynamics (QED) and quantum photonics and has been a subject of great interest for 
superconducting qubits da , atomic [SHS] and quantum dot jl] systems. In addition to elec¬ 
tronic and optical frequencies, strongly coupled light-matter systems can also be obtained 
for quantum well intersubband resonators coupled to mid-infrared [7] and THz [HI [9] frequen¬ 
cies. When a transition is strongly coupled to a cavity resonance, the bare frequencies of the 
uncoupled system shift to new frequencies. The frequency shift depends on the strength of 
the coupling and can be explained in terms of periodic exchange of energy through vacuum 
Rabi oscillations mm or the coupling between the oscillators in a strongly dispersive 
system [12]. Since frequency shift is proportional to light-matter coupling, switchable and 
tunable strong coupling is of interest for hlters and modulators. Over the years, artihcial 
structures of sub-wavelength sizes with novel electromagnetic properties, commonly known 
as metamaterials (MMs) [TH] , are used to conhne and enhance the free-space electromagnetic 
radiation in the near held and hence to study the light-matter interactions in both strong 
and ultra-strong regimes. 

Recently, ultrastrong coupling experiments have been performed, where the cyclotron 
resonance of the two dimensional electron gas (2DEG) is coupled to the THz split ring 
resonators (SRRs) [TH [T5] . complementary SRRs [T6| and superconducting complementary 
metasurfaces HZI by externally tuning the magnetic field. However, electrical tnnability is 
more preferable for extending the potentials of light-matter coupling into realizable devices. 
In 2012, Gabbay et al. ra have modelled the coupling of metamaterials to electrically 
tunable ISRs in a square quantum well for mid-infrared frequencies, later demonstrated 
experimentally by Benz et al. [T9H2T] . Due to local bending of the incident electric held in a 
direction parallel to the growth axis (polarization selection rule) around the metamaterials, 
they can be employed to couple the free-space THz radiation to excite the intersubband 
transitions in parabolic quantum wells [22H21| and modulation doped multi-quantum wells 
[2S|- High electron mobility transistor (HEMT) structures, on the other hand, form the 
building blocks of most modern high-speed electronic circuits. Monolithic integration of 
MMs with HEMT permits amplitude modulation of THz radiation at MM resonances up 
to few MHz [26] in the linear light-matter regime. Such integrated devices not only show 
potential application towards electrically tunable THz devices, but also in the ultrastrong 
light-matter interaction regime. 


2 


In this contribution, we exploit the wide tuning possibility of ISR in a single triangular 
quantum well by an external bias in contrast to the square well, by driving the device from an 
uncoupled to a strongly coupled regime. We present experimental results on strong coupling 
of the 2DEG ISR in a modulation-doped HEMT structure to the fixed cavity resonance of a 
THz MM by electrically tuning the ISR with a high quality epitaxial, complementary-doped 
and transparent electrostatic gate The tuning mechanism is attributed to the quantum- 
confined Stark effect |2H]- The use of complementary epitaxial gates is advantageous since 
the low vacuum band-offset leads to greater electrical tunability. 

Designing the Metamaterial-HEMT Device 

Intersubband transitions in semiconductor quantum well can be designed to cover the wide 
infrared region of the electromagnetic spectrum [221EO] • The device is based on intersubband 
transition in 2DEG, formed in a GaAs/ARGai-ajAs heterojunction of a HEMT structure, 
fig-B a). The HEMT consists of a heterostructure that is delta-doped in the ARGai_a,As 
layer, 50 nm from the heterojunction. The electrons from the delta-doped donors diffuse to 
the GaAs layer. This results in the presence of a strong built-in space charge fields between 
the carriers in the GaAs and donors in the ARGai-^jAs. The potential near the hetero¬ 
junction can be approximated as a triangular potential with a constant slope in GaAs layer 
from the space-charge field and a vertical barrier from the GaAs/ARGai-ajAs conduction 
band-offset. The electronic motion perpendicular to the surface is quantized, resulting in 
the formation of quasi-two-dimensional-subbands. By varying the external bias, the slope 
of the triangular potential along with the transitions energies can be tuned. In contrast to a 
square with abrupt barriers on either side, a triangular well with gradual slope (on the GaAs 
side) will have a much greater effect on tuning the subband spacings by electrical bias. The 
external bias is applied by making contacts to the 2DEG and a heavily doped p-type region 
on the surface which forms an epitaxial gate. The complementary doping of the epitaxial 
gate leads to a low vacuum band-offset (typically 0.5 eV), thus a lower depletion of carriers 
and hence a wide voltage tunability of the device (or intersubband levels). The thickness 
of the epitaxial gate is only 20nm which is much less than the skin depth. This makes the 
epitaxial gate transparent to THz radiation. Further details on the layer sequence of the 
sample are described in the methods section. 

Arrays of double-sided Au split-ring resonators (SRR) are deposited on top of the de- 
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vice, shown in fig. [^b) and (c). The SRRs have sub-wavelength sizes and are designed to 
have a cavity resonance comparable to the first intersubband transition in the 2DEG. The 
resonators are similar to electrical-inductor-capacitor (ELC) circuit, where the inductive re¬ 
gions correspond to the loops of the SRRs and the capacitive regions correspond to the gaps 
in the loops of the SRRs. As shown in fig. id). the double-sided SRRs consists of two single 
SRRs that are mirror images of each other for reflections along the central metal line H 
The current travels in a clockwise direction for one of the single SRRs and in anticlockwise 
for the other. Thus the induced magnetic fields from the two SSRs will average to zero and 
the double SRR will only couple to external electric fields through the capacitive regions. 

Einite difference time domain (FDTD) simulations using a commercial software package 
IHH is used to calculate the electric field distribution and the frequency response of the MM 
modes. The simulations are carried out for a single unit cell with periodic boundary condi¬ 
tion mx — y plane and open boundaries (with spaces) in .2 plane. The electric field direction 
of the incident signal is polarized along y direction. A schematic of the arrangement used for 
electromagnetic simulation along with the dimensions of MM are shown in fig. [^d). The in¬ 
plane field distribution is strongly enhanced over sub-wavelength volumes 

around the capacitive regions of the double-sided SRRs, fig. [^a) and (b). Intersubband 
transitions only couple to electric fields polarized in the growth direction perpendicular to 
the surface. Although the fields between the capacitor sections are predominantly in-plane, 
there is a significant out-of plane field from fringing effects on the surface [32]. This can be 
seed in fig. [^c), where the electric field distribution, {\Ez\) in the growth direction along the 
cut, indicated by the white line in fig. [^a) and (b) is shown. The equivalent electrostatic 
capacitance is an important parameter of the microcavities which play a key role in strong 
light-matter coupling. This is due to the fact that the Rabi frequency is directly propor¬ 
tional to the equivalent capacitance [33|. 

2DEG Intersubband Resonances 

In order to characterize the energy eigenvalues of the 2DEG confined in a triangular potential 
well, a self-consistent Schrodinger-Poisson equation is solved using the ID Poisson solver [3l] . 
The solution is carried out for different bias applied to the gate. The dependence of the 
intersubband spacings of two resonances, Eq Ei and Eq E 2 are plotted as a function of 
the electric field strength in fig. |^a), where Eq, Ei and E 2 corresponds to the ground, first 
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and second subbands respectively. By tuning the gate voltage, the confinement potential and 
hence the electric field strength, F can be altered. Under the triangular well approximation 
[SS] , the shift of the ISR to higher energies is given by [SH] : 



( 1 ) 


where i (0, 1, 2, ...) represent the indices of the subbands, m* is the effective mass of 
electrons in GaAs and e is the electronic charge. The normalized oscillator strengths, fij, 
are calculated from the transition matrix elements, {j \ z |i), corresponding to the ISRs by 
using the following relation: 



( 2 ) 


where |(j| |i)| is given by [36] : 



(3) 



. The quantity, L is defined as the electric length, given 
The normalized oscillator strength for Eq Ei transition is 


calculated to be 0.73 while that of Eq E 2 transition is 0.12. According to the sum rule, 
the normalized oscillator strengths of all transitions sum up to 1, which indicates transition 
to higher levels are very weak. Under the framework of triangular well approximation, the 


oscillator strengths is calculated to be independent of the electric field strength (fig. |3ta)). 
This is unique compared to that of the square potential well as shown by Benz et al. m, 


where there is a small contribution of the quantum-confined Stark shift to the oscillator 
strengths in the presence of external bias. Due to asymmetric confinement, the ISRs can be 
tuned with a magnitude one order higher as compared to an equivalent square well. For a 


triangular well with an effective thickness of 20 nm, if the electric field is tuned by 4.5 x 


10^ V/cm, the ISR {Eq Ei) can be tuned by 35 meV, while for an equivalent square well 
for the same change in electric field the ISR is tuned only by 3 meV. The experimentally 
observed ISRs are blue-shifted due to resonant screening from higher subbands, also known 
as the depolarization effect jSSHST] and is explained in the methods section. 

Results 

Characterization of metamaterials 
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The MMs are characterized by THz time-domain spectroscopy (TDS) at room temperature 
under purged conditions, to avoid the water absorption lines. The experimental details 
of the set-up are explained seperately in the methods section. At room temperature the 
thermal energy, ksT (= 25 meV), is greater than the subband energy spacings. This results 
in the occupation of higher subbands, which consequently prevents us from observing the 
ISR (Eo El) from the 2 DEG. Hence under this condition, the response from the sample is 
purely due to MMs. The transmission is measured for two orientations of the MMs, fig. |^b). 
In the inset, the spectra corresponding to the two orientations of MMs are shown. The elec¬ 
tric field of the THz radiation couples to the MMs only when they are oriented parallel to 
the optical table, marked as MetamateriaLO° in fig. |^b). The normalized transmission is 
plotted by dividing the transmitted intensity with respect to the MM orientation when the 
fields do not couple. Two strong resonances are observed in the transmission spectrum at 
0.34 THz and 1.71 THz. The quality factors of the resonance frequencies at 0.34 THz and 
1.71 THz are found to be 0.58 and 2.01 respectively. These resonances have different origin 
described as follows: The lower frequency at 0.34 THz is due to the LC equivalent circuit 
where the current circulates in the inductive part and the field is concentrated between 
the capacitor plates. The resonance at 1.71 THz is due to the half wave resonance j38] 
arising from the sides of MM. There are additional losses arising from the epitaxial layer, 
damping and change in the effective dielectric constant of the layers with the change in 
temperature. From the quality factor of the resonance at 1.71 THz, we estimated the losses 
to be around 0.85 THz. On cooling the sample, the cavity resonance will be blue shifted 
according to j39]: (Uq = + 7 ^, where coq is the bare cavity resonance without any losses, 

Ur is the resonance of the MMs including losses and 7 is proportional to the ohmic losses 
in the metal and radiative damping in the cavity. The resulting resonance frequency with¬ 
out losses is calculated to be 1.9 THz, which is very close to the observed resonance at 4.2 K. 

Charging spectroscopy of the 2DEG 

The charging of the 2DEG subbands is observed by capacitance-voltage spectroscopy at 
4.2 K, shown in fig. The change in capacitance between the top gate and the bottom 
2 DEG layer is measured as a function of the gate voltage by applying DG -|- AC voltage 
to the gate and measuring the AC component of the current across the ohmics. The steep 
increase of capacitance in the voltage ranging from —0.6 V to —0.4 V indicates that the 
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2DEG layer is filled with charge carriers (green curve). When the sample is illuminated with 
far-infrared (FIR) broad-band source (Hg-arc lamp), the threshold voltage shifts more to 
negative values. This is due to the activation of the donor-exchange (DX) centres [10]. As 
the sample is illuminated longer, more donor atoms from the h-doped AlGaAs layer ionize 
and less bias is required to charge the 2DEG subbands, resulting in the shift of charging 
slope to more negative biases. When the charging spectrum does not shift any more (shown 
by the blue curve), a steady state is reached. The activation of the DX-centres is a result of 
residual near band-gap illumination (after filtering by black polyethylene window) from the 
FIR source. The blue shaded region in fig. marks the voltage range where density-chopped 
infrared transmission measurements are performed. The Schrodinger-Poisson equations are 
solved over the shaded voltage range, corresponding to the electric field strengths plotted 
in fig. I^a). 

Strong coupling of metamaterials with 2DEG 

An FTIR transmission set-up with a rapid scan BRUKER IFS113V interferometer is used 
for the density chopped infrared transmission spectroscopy to observe the strong coupling 
of the MM cavity resonance with the 2DEG ISRs. As observed in the charging spectrum, 
at a bias of —2 V, the conduction band of the 2DEG is pulled above the Fermi level and 
hence completely depleted of carriers. This voltage is used as the reference voltage in the 
chopping-scheme at which sample transmission is recorded. The voltage is then increased to 
— 1.50 V, where only one subband is below the Fermi level and corresponding transmission 
measurements are performed. The voltages are changed alternatively and the respective 
spectra are co-added and averaged over time. This is repeated for several gate voltages, 
keeping the reference voltage same. All the spectra are collectively plotted in a contour-plot 
in fig. I^a). It is observed that at low temperature, the cavity resonance shifted to 2.12 
THz due to lower losses as compared to the room temperature measurements. The Q-factor 
of the MM resonance has dramatically improved to a value of 16.9. At = —1.425 V, 
a clear splitting of the ISR can be observed when the ISR (wqi) of the 2DEG crosses the 
resonance of the MMs at 2.12 THz. The width of the splitting is found to be 0.52 THz. 
Two normalized transmission spectra at = —1.45 V and = —1.40 V, shown by the 
white line in fig. |^a) are plotted in fig. [^b) and (c). The spectra are smoothed in order 
to carry out the spectral deconvolution. Three peaks, marked as (Uoi, ^^meta and u}q 2 can be 
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seen corresponding to the points marked in fig. |^a). As described before, the ISRs are blue 
shifted by 1.2 meV due to the resonant screening from the higher subbands. It is observed 
that at = —2 V, not only the 2DEG channel gets completely depleted of charge carriers 
but also the response from the capacitive element of the MMs is screened due to charges 
in the epitaxial layer. By applying bias on the epitaxial gate below the MM, a significant 
modulation of the intensity (or amplitude) of the cavity resonance is observed [261, HI] • Thus 
in our chopping scheme, the MM resonance at 2.12 THz disappears for biases lower than 
— 1.45 V and begins to appear for biases higher than —1.425 V. In our experiments, we limit 
the bias applied on the gate to —1.3 V. Since with increasing bias the slope of the triangular 
confinement becomes more steeper, the second subband is pulled below the Fermi level, 
which makes the transition scheme different, which is not desirable. 

Discussion 

The 2DEG and MMs in our device can be considered as two oscillators, one of which has a 
fixed frequency and the frequency of the other oscillator (2DEG) is tuned by the gate voltage. 
When the frequencies of both the oscillators are similar, they form a coupled system and 
an anti-crossing phenomenon is observed. Using the common oscillator model, described by 
Gabbay et al. ra. the strength of this coupling is found to be directly proportional to the 
strength of the splitting, described in the methods section. The quantity, fl, defined as the 
coupling strength is given by: 

n=UJ+-UJ. = - Go)]^ (4) 

where Uq is the bare coupling strength, ^meta and ^lo represent the losses in the metamterials 
and ISRs and can be equated to the linewidths of the respective transitions. The coupling 
strength in our experiment is found to be 0.52 THz, while the bare coupling strength is 
calculated to be 0.56 THz. There is no significant change in the two coupling strengths. 
The ratio of the splitting to the sum of full-width at half maximum of both the resonances 
^Aco — «^+Aa)io ~ ^ ~ 1.44) is found to be greater than one. The normalized coupling ratio 
(^) in our experiment is found to be 0.26, which is a clear indication that the observed 
coupling is in the strong regime and is a significantly high value achieved for a single well. 

In conclusion, we show strong coupling of the 2DEG ISR in a heterostructure with a 
triangular confining potential to the cavity resonance of a THz MM by driving the device 




from an uncoupled state to a strongly coupled state via electrical tuning of ISR. Due to 
asymmetric confinement potential, the tuning of the ISR is found to be one order of magni¬ 
tude higher in comparison to an identical square well. The use of epitaxial gate proved to be 
advantageous for the use of the direct MMs and consequently tune the ISR at the same time. 
With a proper chopping scheme, we successfully demonstrated that the measurement can be 
performed in one integrated sample without the need of additional reference sample. This 
is one of the remarkable feature of our device. Moreover, we succeded to achieve ultrastrong 
light-matter interaction by employing a single triangular quantum well in a high electron 
mobility transistor heterostructure, with a normalized coupling ratio of 0.26. 

Methods 
Sample Design 

The sample is grown on a semi-insulating GaAs (100) substrate by molecular beam epitaxy. 
At first 20 periods of AlAs/GaAs (5 nm/5 nm) short period superlattice (SPS) are grown 
to smoothen the surface for growth. A 650 nm thick GaAs layer is grown followed by a 
50 nm Alo. 34 Gao. 66 As spacer layer. This is followed by a Si-5 doping and then a 45 nm 
Alo. 34 Gao. 66 As layer. Another 6 periods of SPS layer comprising of 1 nm GaAs and 3 nm 
AlAs are grown, capped by a 13 nm GaAs layer above which the epitaxial gate |27j is 
grown. The epitaxial gate has a 15 nm thick carbon-doped GaAs layer and 30 periods of 
carbon-5 doped and 0.5 nm carbon-doped GaAs layers, which lead to an effective carrier 
density of = 2 x 10 “^® cm“^. To establish a contact with the 2 DEG layer, 162 nm 
below the surface, the corners of a 5 mm x 6 mm sample are first etched down by 100 
nm (fig. [^b)) and then indium is diffused in a reaction chamber in an inert atmosphere of 
argon and nitrogen. The sample is then mounted on a chip carrier and wires are bonded 
by wedge-bonding. The structures are processed on the epitaxial gate, shown schematically 
in fig. [Jb), by conventional UV-photolithography, evaporation of Gr/Au ( 10/200 nm) and 
lift-off technique. Scanning electron micrographs (SEM) of the MMs and a wire bonding 
are shown in fig. [T](c) and (d). 

THz time domain spectroscopy 

A Ti:Sa laser with 80 fs pulse duration (centre wavelength of 800 nm) and a repetition rate 
of 80 MHz is used to generate THz radiation by exciting an inter-digitated photoconductive 
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antenna processed on a GaAs substrate at an applied DC bias. The experiment is performed 
in a transmission geometry with an NIR-power of 200 mW on the antenna. The detection 
is based on electro-optic sampling of the THz electric field by employing a birefringent, 2 
mm thick ZnTe crystal. The electric field component of the source is in a direction parallel 
to the optical table. Two 90° off-axis parabolic mirrors are used to collimate and focus 
the THz beam on the sample, which is fixed on a rotation mount. Another two parabolic 
mirrors are used to collect the signal from the sample and focus it on the ZnTe crystal. 
The transmission measurements are performed to observe the response of the metamaterials 
alone in the device under purged conditions. 


Oscillator Strengths and Depolarization Shifts 

Under the triangular well approximation, the intersubband energy edges, U*, are directly 
proportional to where F is the electric field strength. Now the oscillator strength for the 
transition between two subbands in a triangular potential well is proportional to the square 
of the matrix elements, which in turn is proportional to the square of the electric length, 
L. The electric length is inversely proportional to This implies that the oscillator 

strength does not depend on the field strength in a triangular well and is a constant given 
by: 


k = u 


TT 


\2/3 12 


2 ^ {U-tk 

-| 2/3 




(5) 


where, ti = — ‘in [i + ^ j2 ' [HE]- For z = 0 and j = 1, /oi = 0.73 and for z = 0 and j = 2, 

/o 2 = 0.12. Taking in to account the allowed transitions to higher subbands {j — 3,4,...), 

the total oscillator strength, J2foj^ sums up to 1. The normalized energy eigenfunctions for 

j 

the ground and the first excited intersubbands under triangular confinement, are given by 


'#>.w = 4rV('7-Av) 


( 6 ) 


VL eFLj ’ 

1 /3 

where. A* are the normalization constants given by Ai' [ti] and L = (h^ j2m*eF'j . In¬ 
cluding the resonant screening effects, the intersubband energies are blue shifted by 


Fio — Fioy (1 -l- CKii), 

where an is the depolarization factor, which according to Ando is given by: 

ine^n2d ^ 1 

ttii —-<^ 117 ;—, 

eeo Fio 


(7) 


( 8 ) 
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where n 2 d is the 2 D carrier density; e is the dielectric constant of GaAs; cq is the permittivity 
of free space. Uujiq is the energy spacing corresponding to the transition, Eq —?• Ei and the 
overlap factor, An has the dimensions of length and is given by [12]: 


OO 


2 


All ~ J j (A) Ao (A) 


(9) 


0 Lo 


which can be written as: 



( 10 ) 


The carrier density is however not a hxed value and changes with the change in bias. Besides 


the bias also changes the slope of the conhnement thus changing the overlap-integral (An). 
At a bias of —1.4 V, corresponding to a carrier density of 2 x 10^^ cm“^, the depolarization 


shift is calculated to be 1.2 meV. The details of the calculation is shown in the methods 
section. The influence of resonant screening becomes important when the level spacings are 
smaller and the carrier densities are high. 

Coupling Strength 

The system is assumed to be comprised of two oscillators which are tuned from uncoupled 
to strongly coupled regime by tuning the ISR of the 2 DEG with an external gate bias, 
making use of quantum-conhned Stark effect. Following the model of coupled oscillators, as 
developed by Gabbay et ah [T 8 |, the coupling between the two resonators can be described 
by a simple 2 x 2 matrix: 



( 11 ) 


Do/2 uio — ! 


where (Uio, are the resonances and ^lo, ^meta are the damping factors corresponding to 

ISR and metamaterials; Dq is the bare coupling strength of the oscillators. The matrix can 
be diagonalized to get the upper and lower eigenvalues of the coupled system as: 



Subtracting the frequencies, we obtain: 



( 13 ) 
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At the anticrossing point, both the resonators start oscillating with equal and real frequen¬ 
cies, the coupling strength, f], can be expressed as: 


^ — UJ- — 2]j-Q q — [{^rneta ~ ^lo)]^ • (14) 

From the two (upper and lower) branches of the transitions, the coupling strength, fl, is 
found to be 0.52 THz. The damping factors corresponding to the resonances are ^rneta = 0.13 
THz and = 0.23 THz. The bare coupling strength, is calculated to be 0.56 THz 
which similar to the coupling strength indicating strong coupling with low losses. 
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FIG. 1. Device design, (a) Conduction band edge, Ec, of our device when a certain bias is applied 
at the gate such that the ground state is below the Fermi energy. The ground wavefunction is also 
plotted within the well, (b) Schematic layer sequence of the device and the corresponding electric 
field direction of the incident infrared radiation which couples to the resonator. The bottom surface 
of the device is wedged at an angle of 3° in order to avoid the Fabry Perot fringes in transmission 
measurements. SEM micrograph showing (c) metamaterials and wire bonding, (d) metamaterial 
unit cell. 
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At z = -0.162 |am 



FIG. 2. FDTD simulations of metamaterials. In-plane electric field distribution, 
(yW -h l^yP), at (a) z = 0 pm (on top of the surface) and (b) z = —0.162 pm (in the 
2DEG layer, 162 nm below the surface), (c) The field distribution, \Ez\^ in the growth direction 
along the cut shown by the white lines in (a) and (b). The black dotted line represents the position 
of the 2DEG. (d) The schematic representation of the metamaterial structure and the electric field 
component of the polarized excitation source for the TDS measurements. All dimensions are in 
the units of pm. The metamaterials are arranged in an unit cell of 166 pm x 182 pm. 
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FIG. 3. ID Poisson simulation and characterization of metamaterials, (a) Dependence of 
the intersubband resonances and the oscillator strengths for two intersubband transitions, |0) ^ |1) 
and |0) ^ |2), on the electric field strength and hence the bias applied on the gate. The oscillator 
strengths do not vary with the gate voltage and are fixed for each transition, which is a characteristic 
feature of the traingular potential well. Inset: Band-schematic of the first three subbands in a 
2DEG with respect to the Fermi level at an intermediate gate voltage where only the ground 
subband is below the Fermi level, (b) The transmission of the incident polarized electric field for 
the horizonzal orientation of the metamaterial (MetamateriaLO°) normalized to the transmission 
for the vertical orientation of the metamaterial (MetamateriaL90°). The green spectra are results 
of deconvolution and the red curve is the reconstructed transmission spectrum. Inset: Spectra 
corresponding to the two orientations of the metamaterials. The schematic shows the orientations 
of the meta-structures with respect to the electric field direction in the THz TDS measurements. 
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FIG. 4. Capacitance-Voltage spectroscopy. Dependence of capacitance on the voltage applied 
to the gate. With illumination, visible part of the source saturates the DX centres and hence the 
threshold voltage shifts to more negative values. When a quasi-stable condition is reached, and the 
thereshold does not shift any further, all donors are saturated. The blue shaded region highlights 
the range of voltage over which the FTIR transmission spectra are taken. 


19 







3.75 

3.50 

3.25 

'n' 

f 3.00 


2.75 
2.50 
2.25 
2.00 

1.75 



-1.500 -1.475 -1.450 -1.425 -1.400 -1.375 -1.350 


Gate Voltage (V) 



Frequency (THz) 


FIG. 5. Observation of strong coupling, (a) Strong coupling is observed in the transmission 
spectra by tuning the gate voltages for the | 0 ) ^ | 1 ) transition with the metamaterial resonance 
at around 2 THz, while no coupling is observed for the other intersubband transition, |0) ^ |2). 
The coupling of the cavity resonance with the 2DEG intersubband resonance is seen as a split at 
the cross-over point between the two resonances with a split-gap of 0.52 THz. The normalized 
transmisison spectra at two voltage points, (i) = —1.45 V and (ii) = —1.40 V are shown 

in (b) and (c) respectively. Each spectral deconvolution shows three dips (cjqi, ojmeta and 6 ^ 02 ) 
corresponding to the points shown in (a). The blue dashed lines show the shift of the respective 
transitions. 
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